Asteroid 2007 WD5 will not impact Mars on January 30! 
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ABSTRACT 

The Monte Carlo method of the nominal orbit clonning was applied to the case of 2007 WD5, 
the asteroid from the Apollo group. Calculations based on 33 observations from the time interval of 
2007 11 08 - 2008 01 02 showed that the asteroid will pass near planet Mars at the minimum distance of 
10.9 ±2.9 Rjviars) what implies that probability that 2007 WD5 strike the planet decreased to the value 
of 0.03% from the value of about 3-4% previously announced by NASA. The additional observations 
taken on January 3-9 reduce further the asteroid's impact chances, effectively to nil: the asteroid will 
pass near planet Mars at the minimum distance of 8.4 ± 1.1 Rjviars- 

1. Introduction 

The asteroid, known as 2007 WD5, was discovered on November 20, 2007 by Boattini of 
the NASA-funded Catalina Sky Survey on Mount Lemmon, near Tucson (Arizona, USA) 
using a 1.5 m telescope. The size of asteroid was estimated as ~50 m. 

Next day, Chesley and Chodas (2007) from the NASA/JPL Near-Earth Object Program 
Office announced that newly discovered asteroid which passed close to the Earth in November, 
will pass very close to Mars in late January, and there is a chance that it could hit that planet. 

On December 28, 2007, Yeomans et al. (2007) informed on Web page that three pre- 
discovery observations of asteroid, taken on November 8, 2007 were located by Puckett in the 
archive of the Sloan Digital Sky Survey II at the Apache Point Observatory. They refined 
the orbit of 2007 WD5 and found that the probability of Mars impact on January 30, 2008 
increased to 4%, whereas the uncertainty region of the asteroid position during the Mars 
encounter has shrunk to 400 000 km along a very narrow ellipsoid. They stressed that the 
uncertainty region of the asteroid position during the Mars encounter could shrink more - to 
the region which will no longer intersect with the planet surface. 

The last (at the moment of writing this text) report from NASA/JPL (Chesley et al. 
2008) - based on the additional observations of the asteroid between Dec. 29 and Jan. 9 - 
significantly diminished the asteroid chances to strike Mars to approximately 0.01%. Chesley 
et al. (2008) estimated that asteroid 2007 WD5 most likely will pass about 8 Mars radii, 
TMars =8RMars5 from the planet's center on Jan. 30 and with 99.7% confidence, the asteroid 
should be no closer than 1.2RMars from the planet's surface. 

According to NeoDys (2008) calculations based on the observation arc from 20071108 
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to 2008 0109 the asteroid will pass near Mars at the nominal distance of 2.204 • 10"'^ AU 
(9.7RMars)) whereas the minimum possible distance is 8.815 • 10~^AU (3.9RMars)- These 
estimates suggest that asteroid 2007 WD5 will pass Mars even in the greater distance than 
was suggested by Chesley et al. (2008). 

The aim of this note is to present how our method of the nominal orbit clonning allows 
to determine the set of impact orbits as well as the probability of this impact for the case of 
2007 WD5. 

2. The method 

The present investigations are based on the archive positional observations taken from 
the NeoDys (Near Earth Object - Dynamic Site, University of Pisa, Italy) publicly available 
on the Web at 

|http: / / newton.dm.unipi.it / neodys / mpcobs / 2007WD5.rwo[ The whole observational mate- 
rial contains 46 observations covering the time period from November 8, 2007 to January 9, 
2008. 

The equations of asteroids motion have been integrated numerically using the recurrent 
power series method (Sitarski (1989), (2002)), taking into account the perturbations by all 
the nine planets and by the Moon. The planetary coordinates were taken from the Warsaw 
numerical ephemeris DE405 / WAW of the Solar System, consistent with a high accuracy with 
the JPL ephemeris DE405 (Sitarski 2002). 

We derived the nominal orbits using the least squares orbit determination method for four 
different sets of residuals given in the second column of Table [T] The respective solutions are 
denoted as A, B, C, and D, where solutions B and C differ only in a number of residuals taken 
into account for the orbit fitting: three residuals of significantly larger values in comparison 
to the rest of data were excluded in the case C. The nominal rms's are given in column 4 of 
Table [H 

The sample of 10 000 clones of each nominal orbit was selected randomly according to the 
normal distribution in the 6D-space of orbital elements (the method is described in details 
by Sitarski (1998)). In the column 5 of Table [l] upper limits for the rms of the true (but 
unknown) orbit are given for the confidence level of 99%. 

The orbit clones from the samples A-D were integrated forwards up to the close encounter 
with Mars on January 30. 

2. Results 

We obtained 220 impact orbits for the Mars encounter on January 30, 2008 by numerical 
integrations of 10 000 orbital clones in case A. The distribution of minimum asteroid distance 
from Mars are presented for the solution A by black histogram in Fig. [T| Thus, the impact 
probability of 2.2 • 10^^ was calculated for this arc of observations. However, the refinement 
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Table 1: Description of four orbital solutions for 2007 WD5 and the derived minimum distance 
for asteroid close encounter with Mars on January 30, 2008. To estimate the value of rmstme 
(column 5) the confidence level of 99% was assumed. 



Solution 


Observational 


Number 


Mean 


rmstruc 


Minimum 


Impact 




interval 


of 


res. 


not greater 


distance 


probability 






res. 




than 


[R-Mars] 


[%] 


A 


20071108 - 200712 31 


64 


0'.'521 


0'.'585 


11.9 ±5.9 


2.2 


B 


20071108 - 2008 0102 


68 


0"512 


0'.'572 


13.0 ±4.6 


0.4 


C 


20071108 - 2008 0102 


65 


0"308 


0'.'346 


10.9 ±2.9 


- 0.03 


D 


20071120 - 2008 0109 


89 


0'.'268 


0'.'293 


8.4± 1.1 


- 0.0 



of the asteroid's orbit using also the observations taken on January 2 significantly reduces 
the chance of impact on Mars. Cyan, blue and magenta histograms in Fig. [T] represent the 
distribution of minimum asteroid distance from Mars for solutions B, C and D, respectively. 
One can see, that for the solution C the impact probability decreases to the value of 3 • 10~^ 
(we found three impact orbits among the sample of 10 000 clones), and for the case D - the 
impact seems to be impossible. 

Gaussian distribution of rMars and the impact probability are given in the columns 6 and 
7 of Table 1 and in Fig. [T] Projections of the 6D ellipsoid of orbital elements onto the a-e 
plane for the solutions A-D are shown by red points in Fig. [2]- [3} The shrinking of the 
uncertainty ellipsoid with the refinement of 2007 WD5 nominal orbit and its moving away 
from the impact region are clearly visible. For the longest arc of observations (case D) the 
asteroid will pass near planet Mars at the minimum distance of 8.4 ± 1.1 RMars (magenta 
distribution in Fig [T]) what implies that with 99.7% confidence, the asteroid should be no 
closer than 4.1 RMars from the planet's surface. 

Thus, our calculations indicate that the possibility of asteroid 2007 WD5 collision with 
Mars seems to be ruled out. 
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Figure 1: Distribution of the minimum distance of the asteroid 2007 WD5 from the centre of 
Mars in 2008 01 30.5 derived for the samples of 10 000 clonned orbits. Black, cyan, blue and 
magenta distributions correspond to the solutions A, B, C and D, respectively. 
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Figure 2: Projection onto the a-e plane in the 6-dimensional space of possible osculating 
orbits of 2007 WD5 obtained for solution A (upper panel) and solution B (lower panel) Sample 
of 10 000 clonned orbits is given by red points, while the impact orbits are given as black 
crosses. The plot is centered on the values of semimajor axis, ao, and eccentricity, eo, of the 
respective nominal orbits (epoch: 2007 10 27). 
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Figure 3: Projection onto the a-e plane in the 6-dimensional space of possible osculating 
orbits of 2007 WD5 obtained for solution C (upper panel) and solution D (lower panel). 
Sample of 10 000 clonned orbits is given by red points, while the impact orbits are given as 
black crosses. The plots are centered on the values of semimajor axis, ao, and eccentricity, 
Co, of the respective nominal orbits (epoch: 2007 10 27). 
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